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Abstract 

The two-matrix model is defined on pairs of Hermitian matrices (Mi, M2) of size n x n 
by the probability measure 

^- exp (Tr(-y(Mi) - W(M 2 ) + rMiM 2 )) dM 1 dM 2 , 

where V and W are given potential functions and r G R. We study averages of products 
and ratios of characteristic polynomials in the two-matrix model, where both matrices Mi 
and M2 may appear in a combined way in both numerator and denominator. We obtain 
determinantal expressions for such averages. The determinants are constructed from several 
building blocks: the biorthogonal polynomials p n (x) and q n (y) associated to the two- matrix 
model; certain transformed functions P n (w) and Q n (v); and finally Cauchy-type transforms 
of the four Eynard-Mehta kernels Ki,i, -Ki,2, ^2,1 and K 2 ,2- In this way we generalize known 
results for the 1-matrix model. Our results also imply a new proof of the Eynard-Mehta 
theorem for correlation functions in the two-matrix model, and they lead to a generating 
function for averages of products of traces. 

Keywords: Two-matrix model, average characteristic polynomial, Eynard-Mehta the- 
orem, biorthogonal polynomial, multiple orthogonal polynomial, Riemann-Hilbert problem, 
determinantal point process. 

1 Introduction and statement of results 
1.1 The two-matrix model 

Let 7i n be the space of Hermitian n by n matrices. The two-matrix model is defined on ~H n x % n 
by the probability measure 

exp (Tr(-y(M x ) - W(M 2 ) + T M 1 M 2 )) dM x dM 2 , {M 1 ,M 2 ) e U n x H n , (1.1) 

where V and W are given real- valued potential junctions, where r E K is the coupling constant, 
Z n is a normalization constant, Tr denotes the trace and dM\dM 2 is the Lebesgue measure on 
% n x Hn, see [53]. The matrix model (|l.ip is probably the most popular instance of a two-matrix 
model in the literature, see e.g. [IJ|![l0l[IlJ[I3[21J[2![23[2^ among many 

others. We may mention however that other kinds of two-matrix models exist as well, see e.g. 
IT21 and the references therein. 
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The normalization constant Z n in (jl.ll) is called the partition function. It is defined as the 
total integral 

Z n = / / exp (Tr(-F(M 1 ) - W(M 2 ) +tM 1 M 2 )) dM x dM 2 . (1.2) 

We assume that the potential functions V and have sufficient increase at infinity in order 
that the integral (jl.2D is convergent. Throughout this paper, we will assume that V and W are 
polynomials of even degree with positive leading coefficients. 

Due to the Harish- Chandra formula (or Itzykson-Zuber formula) [33], integrals over the two- 
matrix model (jl.ll) can often be rewritten in terms of the eigenvalues of Mi and M 2 . We need 
the following version of this result. Assume that the functions /, g : H n — > C depend only on 
the eigenvalues (A^)" =1 of Mi and (yt*i)" = i of M 2 , respectively, (taking into account eigenvalue 
multiplicities), in the sense that 

/(Mi) = /(Ai,...,A„), g(M 2 ) =g(fii,...,fin), 

where /, g : R" — > C arc symmetric functions of n variables. Then 

— / / /(M 1 ) ff (M 2 )e Tr ^ y ( Ml )" w/ ( M2 ) +TMlM2 ) dhh dM 2 

4T-T /(Ai ' • • ■ ' ■ • ' " n) n (e- V ^e- w ^A 

Zj n J —OO J —OO i—\ 

n 

x A(A)A( M ) det(e TA ^ )£ =1 J[ (dX, d^) , (1.3) 

i=l 

for a new normalization constant see [53]. The factors A (A) and A (/it) in (ll.3|) denote the 
Vandermonde determinants Yii<i<j< n (^j ~ ^») an< ^ llKi^^fe ^Z 1 *) respectively. In (|1.3[) we 
also assume that / and g are such that the integrals converge. 



1.2 Biorthogonal polynomials and kernels 

It is known that the two-matrix model (|l.lj) can be analyzed by means of biorthogonal polyno- 
mials and kernel functions that we describe below [23] [5H [33] . We will also introduce certain 
(non-standard) Cauchy-type transforms of these objects, which will be needed in the statement 
of our main theorems. 

The monic biorthogonal polynomials {pi(x)}°^L and {Qj{y)}jLo are defined as the monic 
polynomials 

Pi (x) = x i + o(x i - 1 ), qj (y) = y j + o(y j ~ 1 ), 

that satisfy the orthogonality relations 

/ Pi (x) qj (y)e- V ^- W ^+ Tx y dxdy = 0, i^j. (1.4) 

-co J — oo 

These polynomials exist and are unique, see [23) . 

We define the transformed functions Pi and Qj by 

/oo 
Pi (x)e- V ^- W ^ +Tx y dx, (1.5) 
-oo 

/oo 
^jeWdfH^ (1.6) 
-oo 
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and we let hf be defined by 

h\ = I I Pl (x) qi {y)e- V ^- W ^+ Tx y dx dy. (1.7) 



OO rOQ 



oo J — OO 



The integral fjl.7p is positive, see [23] . 

The Eynard-Mehta kernels #1,2) #2,1 and #2,2 are defined as pMl [T5] 

n—l ^ 

^1,1(^17^2) = j^pi{xi)Qi(x2), (1.8) 

i=0 * 
n— 1 ^ 

#1,2(0:, 2/) = ^2j^p l (x)q l (y), (1.9) 

X"a,i(l/,x) = h£ ^P^Q^ - e -v(*)-w(»)+rx« } (1 . 10 ) 

#2,2(1/1,2/2) = ^T2^(yi)ft(y 2 ). (1.11) 



2=0 1 

It is known that the correlation functions for the two-matrix model (jl.ip can be written as deter- 
minants of matrices built out of the four Eynard-Mehta kernels; see [23] and see also Section l3~21 
below. 

In what follows, we will need certain Cauchy-type transforms of the above objects. First we 
define the Cauchy transforms Pi(w) and Qj(v) of the functions Pj and Qj in (ll.5p ~ p.6p : 

p H = r2W dv= r rM^-mn^r, 

Qj<!°) = / f3^= / qAv)e dZdr). (1.13) 

J— oo ^ S J— ooJ — oo ^ S 

We will also need the following Cauchy-type transforms of the Eynard-Mehta kernels: 

1 f°° x - £ 

K 1A (x,v) = / ^K 1A (x,0dC, (1.14) 

K h2 (x,y) = K 1<2 (x,y), (1.15) 
£ a ,iK») = /°° f p K24 w ,0 n dgd-, (1.16) 

K 2 , 2 (w,y) = — -— / ^-K 2t2 (r),y) dr,. (1.17) 

See Section T2.4I for several properties of the kernels i,j = 1,2. Note that in the above 

formulas, we consistently use and '77' to denote integration variables which play the role of an 
'x- variable' or a 'y-variable' respectively. 

The following lemma gives some alternative expressions for the kernels #i,j, i,j = 1,2. 

Lemma 1.1. (Summation formulas for the kernels Kij:) We have the following equivalent 
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expressions for the kernels Kij in (|1.14[) — (jTTTTJ) : 



n-1 



K ly2 {x,y) = ^2^ I p l (x)q l (y), (1.19) 
K 2A (w,v) = EpflWQiW - / p— - TT— ^ gdi?, (1-20) 



v i=0 



ft? / 7-oo7-oo (w-»?)(w-0 

^ 2)2 Ky) = fr'p^y) (i.2i) 



hf 



w-y 



where we use the notations in (|1.12p - (|1.13[) . 
Lemma 11.11 is proved in Section 12.11 

1.3 Average characteristic polynomials 

Let the matrices Mi,M 2 G "H n be distributed according to the two-matrix model (jl.ip . The 
average characteristic polynomial of the matrix Mi is defined by (the notation in the left hand 
side is explained below) 

pHACO]^ . = 1 f r det{x j n _ Mi) ^Ti( — V (Mi) — W(M 2 )+tM 1 M 2 ) dMi m ^ (L22) 

where /„ denotes the identity matrix of size n. Similarly, the average characteristic polynomial 
of the matrix M 2 is defined by 

p[o,i,o,o]( y) := — ( f det(yl n - M 2 ) e ^(-v(M 1 )-W(u 2 )+rM 1 M 2 ) dMi dM% (1 23) 

In this paper, we consider more general versions of (|1.22l) (|1.23p . We are interested in averages 
of products and ratios of characteristic polynomials, where both matrices M\ and M 2 may appear 
in a combined way in both numerator and denominator. Our object of study is of the form 

p[i,j,K,L] . . . ( Xj ; j/ x , . . . , j/j; ui, . . . , ujj-; tWi, . . . , tWi) := 
J_ f [ nj=i det(xA - Mi) n/ = i detain - M 2 ) eTr( _ v(tfi) _ ff(M2)+TMiM2) ^ ^ 
./«„./«„ nf=i det(u fc / n - Mi) nf = i det(^J„ - M 2 ) 

(1.24) 

Here /, J,K,L G NU {0} are nonnegative integers and we assume that Xi, . . . , Xt, J/i, . . . , J/j G C, 
that vi, . . . , Wi, . . . , wl G C \ R, that the numbers in the set (x\, . . . , x/, Vi, . . ■ , vk) are 
pairwise distinct, and similarly the numbers in the set (y\, . . . ,yj, w±,...,wl) are pairwise 
distinct. The numbers Xi,yj,Vk,Wi are sometimes called external sources |S]. 
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Due to (|1.3p . the expression (|1.24l) can be written alternatively as 
p [i,j,k,l] (x u . . . )Xl ;yi, . . . ,yj;vi, . . . ,v K ;tvi, • ■ ■ ,w L ) 

1 r r -q I nf=i o* - A ») n/=i (% - 



71 7* 

[] (e- v ^e- w ^) A(A)A(/i)det(e^^')L-=i II ( dA * ■ ( L25 ) 



x 



In principle, an expression for (jl.24l) can be obtained as follows: if we define the perturbed 
potentials 

V(x) = V(x) - log ( Uj=l^A ) , W(y) = W(y) log ( H ; ^ h , 
then (jl.24p is just the scaled partition function corresponding to these perturbed potentials: 

P\*' J,K ' L] {xi, . . . ,xr,yi, ■ ■ ■ ,yj',vi, . . . ,v K ;wi, . . . ,w L ) = (1.26) 

with 

Z n := J J exp(Tr(-V{M 1 )-W(M 2 )+TM l M 2 fj dM Y dM 2 . 

Our interest, however, is to express (jl.24l) in terms of the unperturbed potentials V and W. In 
this way we will obtain analogues of known results for the 1-matrix model, which go back to 
Brezin-Hikami Q2], Fyodorov-Strahov [HUES] and others [61 151 [TB ] 125 1 ■ 

To start with, we consider the case where only one external source is present. 

Theorem 1.2. (Average characteristic polynomials of M\ and M 2 :) For the average character- 
istic polynomials in (|1.22p ~ (|1.23p . we have 

P Wfi] {x)=Pn[x)i 
Pt im (y) = qn{y), 

where p n and q n denote the biorthogonal polynomials in (|1.4[) . 

Theorem 1.3. (Average inverse characteristic polynomials of M\ and M 2 :) With (I,J,K,L) 
equal to (0, 0, 1, 0) or (0, 0, 0, 1) in (fl~2l| . we have 



Pj [o,o,i,o] (u) = 



h 2 

n n-l 



pP.o,o,i] (t|;) = ^H j 



where we use the notations in (|1.7I) and (|1.12l) - (jl.l3|) . 



Theorems 11.21 and 11.31 will be proved in Section 12.21 We note that the proofs will be quite 
different from the proofs of the corresponding statements for the Hermitian 1-matrix model, see 
[5J [TBJ [T7J [2H1 • The latter proofs typically use a Heine-type representation of the orthogonal 
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polynomials on the real line in terms of a moment (Hankel) determinant [37) . Such formulas 
also exists in the biorthogonal case [33] , but it seems that a straightforward generalization of the 
arguments for the 1-matrix case is not available. Instead, our proofs of Theorems 11.21 and 11.31 
will rely on the multiple orthogonality properties of p n and q n due to Kuijlaars-McLaughlin (32) . 
Next we consider the case where there are exactly two external sources in f| 1 . 24|) . 

Theorem 1.4. (Two external sources:) With (I,J,K,L) equal to (1,0,1,0), (1, 1,0,0), (0,0, 1,1) 
or (0, 1, 0, 1) in ([TT21]) . we have 

Pg'O'Wfav) = (x-v)K 1A (x,v), 
P^'Wfay) = h 2 n K 1<2 (x,y), 

P^ l \y,w) = (y~w)K 2 , 2 (w,y), 

where we use the notations in (11.7)1 and (|1.18l) — (|1 .21[) , but with the index n in (|1 . 19[) and (|1.20)1 

replaced by n + 1 and n — 1 respectively. 

Theorem ll.4l is a consequence of the more general Theorem II .61 below. 

Finally we consider the general case in (|1.24|) . We first state the result when all external 
sources are in the numerator. 

Theorem 1.5. (External sources in numerator only:) Assume that K = L = in (|1.24[) . i.e., 
(|1.24j) has external sources in the numerator only. Also suppose that I > J. Then 

n J ~ 1 h 2 



pV' J <°'°Uxi xr-vi vt) = lli=0 '" n+t 

r n [Xi,...,xi,y 1 ,...,y J ) - X )U (v-v) 

iii<i<j<i\ x 3 x i) iii<i<j<j\yj y-i) 
(K 1<2 (xi,yx) ... Kx <2 (xx,yj) p n+J (xi) ... p n +i-i(xi)\ 
x det : : : : | . (1.27) 

\Ki t2 (xi,yi) ... Ki t2 (xi,yj) p n+J (x r ) ... p n+J _i(xj)/ 



Here K\_ 2 is as defined in ()1.19|) but with n replaced by n + p where p £ {J, ...,/} may be any 
fixed number. 

Theorem 1 1.5 1 will be proved in Section [2T51 The proof will be based on formulas of Christoffcl 
type [3 |SZ] • A close variant to Theorem 11.51 was obtained by Akemann et al. in [3] . 

In the statement of Theorem II .51 we assume that the kernel K\ 2 is defined in (|1.19D but with 
n replaced by n + p for an arbitrary but fixed number p G { J, . . . , /}. The fact that all values of 
p € { J, ■ . ■ , 7} give the same determinant can be seen by applying elementary column operations 
to (|1.27l) . More precisely, to the jth column of (|1.27l) . j £ {1, . . . , J}, we may add a suitable 
linear combination of the last I — J columns. 

Now we turn to the general case. 



G 



Theorem 1.6. (The general case:) Suppose that I — K > J — L > —n. Then 

Yli,i(V]- w i) 



p[i,j,K,t\ (x x , . . . ,xr,yi, , . . ,yy,vx, ■ . ■ ,v K ;wi, . . . ,w L ) 
Xlidxi- v k ) 



Ui<i<j<i(.^j - x i)Ui< k <i<K( v k-vi) Ui<i<j<Av3 ' -yi)Ui< k <i<L( w k - Wl) 

I ifi,i(a;i,«i) ... if 1,1(2:7, vx) #2,1(^1, «i) ... #2,1(^7, «i) \ 



x det 



#1,1(2:1, «k) 
#1,2(2:1, Vi) 

Ki, 2 (xi,yj) 

Pn+J-L&x) 



#1,1(2:7, vk) 

#1,2(2:7, yi) 
#1,2(2:7,2/./) 

Pn+.J-h{xi) 



#2,1(^1, «k) 

#2,2(wi,j/l) 

#2,2(wi,yj) 
Pn+J-L(wi) 



K2,i( w l,vk) 

#2,2 K, J/1 ) 

K 2 M w L,yj) 

Pn+J-L(W L ) 



\p 



n+I-K-X 



with 



C 



Pn+I-K-l(xi) P n+ I- K -l(wi) 



P 



, (1.28) 



n+I-K- 



-iK)/ 



(-i)(wn£o L "^ 2 



(-i) (W nri, 



t=0 '"n+i' 
i=J-L "n+i! 



ifJ-L>0, 
ifJ-L<0. 



Here the kernels Ki^ are as defined in (|1.18[1 - (|1.21[) &ui im£/i n replaced by n + p where p G 
{J — L, . . . , I — K} may be any fixed number. 
Similarly, if J — L > I — K > —n then 



p [i,j,k,l]( Xi ^ _ _ ,xr,yi, ■ ■ ■ ,yj;vx, . . . , v K ; wi, . . . , w L ) 

YiidXi-Vk) 



c 



UjAvj - w i) 



Hx<i<j<i( x i - x i)Hx<k<i<A v k -m) Hx<i<j<Ayj -yi)Hi<k<i<L( w k - m) 

( #1,1(1!, vi) ... Kx.\{xx,v K ) #1,2(2:1,2/1) ■■• #1,2(2:1,2/7) \ 



x det 



#1,1(2:7, Vx) 
#2,i(wi,«i) 

#2,l(w L ,Wl) 

Q,i+i-k(vx) 



Ki,x{xi,v k ) 
K 2t i(wx,v K ) 

K 2 .i(w L ,v K ) 

Qn+I-I<(VK) 



#1,2(2:7,2/1) 
#2,2(^1,2/1) 

#2,2(^7,2/1) 

q n +i~K{yx) 



\Qn+J-L-l(vx) ... Q n+ 7_7_l(i;K) <7n+7-7-l (?/l) 



#l,2(2:7,2/j) 
#2,2(^1,2/7) 

#2,2(^7,2/7) 

q n +i-K{yj) 
q n +j-L-i(yj)J 



(1.29) 



with 



C = 



{-ifwtfr«-i h l +i , ifi-K>o, 

n+ i, ifI-K<0. 



li=0 "'n+i 



Here #j are as m (|1.18p - (ll.21[) ?rai wi£/i n replaced by n + p for any arbitrary but fixed number 

P e{l-K,...,J-L}. 
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Note that it would be more natural to write the matrix in (|1.28l) in its transposed form; we 
have written it in its present form because of typographical reasons. 

Theorem 11.61 will be proved in Section 12.51 The proof will use Theorem 11.51 together with a 
mechanism to transform an external source in the matrix M\ (or M<x) in the numerator ; into an 
external source in the matrix Mi (or M\ respectively) in the denominator. 

In (11.28)) we assume that the kernels JCy are defined in (I1.18l) - (|1.2ip but with n replaced 
by n + p for an arbitrary but fixed number p G {J — L, . . . , I — K}. The fact that all values 
of p G {J — L, ... ,1 — K} give the same determinant can be seen by applying elementary row 
operations to (|1.28p . More precisely, to the j'th row of (|1.28p . j G {1, . . . , J + K}, we may add 
a suitable linear combination of rows J + K + 1, . . . , I + L, taking into account (|1.18p - (|1.2ip . A 
similar remark applies to (jl.29p . 

In the special case where J = L = in (ll.24p . i.e., when there are external sources in the 
matrix Mi only, then Theorem 1 1 . 61 could be obtained from the Ricmann-Hilbert characterization 
in Section 12.31 together with the results in [TO] ; see also Remark 12.51 in Section 12.31 A similar 
remark holds of course if J = K = 0, i.e., when there are external sources in the matrix M2 only. 
We will not use this approach however in the proofs. 

In Theorem 11.61 we assume that 

min(7 — K,J — L) > —n. (1.30) 

We believe that this assumption is sufficiently mild to cover most of the applications. We also 
expect that similar determinantal formulae may exist if the assumption (jl.30p fails. See [5] for 
the statement of such formulae in a similar context, requiring a lot of case distinctions to be 
made. We will not address this issue in this paper. 

We wish to point out that the above formulas have a lot of similarity with those obtained 
in the physical papers [!7, [8] by Bergere in the context of the normal 1-matrix model (or more 
general complex matrix models), see also 4, 5 . Recall that the normal 1-matrix model is defined 
by the probability distribution 

1 e -Tr(V(M,M*)) dMj (L31) 

on the space M n of all normal n x n matrices M. Here the superscript * denotes the Hermitian 
conjugate, and the potential function V in (11.311) is such that its trace is real-valued and depends 
only on the eigenvalues of M. In this context, one has the following analogue of (|1.3I) : 



— / /(A/) 5 (M*) e - Tr ^( M ' M *)) dM 



Z 



1 />oo POO n n 

= ^T ■ ■ /(Ai, • ■ • , A„)g(AT, ■ • ■ X) J] (e- v ^) A(A)A(A) J] dK, (1.32) 

for functions / and g which depend only on the eigenvalues and which are such that the integrals 
converge. Here the bar denotes complex conjugation. 

In the physical papers [3 [8] , Bergere obtains determinantal formulae for averages of products 
and ratios of characteristic polynomials over the normal matrix model (|1.31|) (and also more 
general complex matrix models). His formulas are written as determinants constructed from 
(bi)orthogonal polynomials and kernels in a similar way as in our results. We should stress 
however that the two-matrix model (11.11) and the normal 1-matrix model (jl.3ip are truly different 
models, exhibiting very different properties. Viewed from this perspective, the similarity of our 
results with those of Bergere [5] may appear rather remarkable. It would be interesting to 
obtain a better understanding on this issue. 
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1.4 Organization of this paper 

The rest of this paper is organized as follows. In Section [5] we prove the main theorems. In Sec- 
tion [3] we give two applications: a proof of the Eynard-Mehta theorem for correlation functions, 
and a formula for averages of products of traces. Finally, Section 0] contains some concluding 
remarks on the more general model of matrices coupled in a chain. 

2 Proofs 

In this section we prove the main theorems. This section is organized as follows. In Section [2.11 
we prove Lemma [PI Section |2"T21 contains the proof of Theorems ll.2l and ll.3l In Section |2~31 we 
prove Theorem 1 1.5 1 In Section l2"H1 we discuss several properties of the kernels Ki.j. In Section l2~51 
we use these properties to prove the general result in Theorem ll.6l 

2.1 Proof of Lemma 11.11 

Equations (11 . 19[) and (|1.20[) are immediate by the definitions. Now let us prove (|1.18|) . Using 
the decomposition x — £ = (x — v) + (v — £) in the integrand in (jl.141) . we obtain 

K 1A (x,v) = r —K 1A (x,0 d£ - -J— K 1A (x,0 d£ 

1 ;K ltl (x,Z) d£-—, (2.1) 



v — £ v — x 



where we used that 



/oo />oc n—1 

K hl (x,Od£ := / ]T p^Q^) d£ 
-co J—oo "i 

/oo p oc n—1 ^ 
/ E T2Pi(^Me- V ^- W(r,)+T ^ # dr, 
-ooJ-oc l=Q H 

n ~ 1 / pOO p OO 

j_Q \J — ooJ — oo 

= 1, (2.2) 

due to the biorthogonality relations and the fact that po(£) = 1< Now by inserting (jTTHJ) in (|2.ip . 
we find that 

x 1 , 1 (x,«) = Et2^w / —7Qi(0de 

■ n "i \J — oo V S J V X 



which is (|1.18p . The proof of (|1.21j) is similar, this time using the following analogue of (|2.1[) : 

/•°° 1 1 

K 2t2 (w,y)= K 2 . 2 ( V ,y) d V . (2.3) 

J.oo io - 77 w-y 

□ 
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Remark 2.1. (Reproducing properties of \ and K% t i\) The argument used to obtain (|2.2[) in 
the above proof shows more generally that 



P (OK hl (x,OdC= P (x) (2.4) 

J — OO 

and similarly 

q{v) K 2Av, y) d v = q{y) (2-5) 



whenever p (or q) is equal to one of the biorthogonal polynomials pi (or qi respectively) with 
i = 0, . . . , n — 1. By linearity, (|2.4[) - ([2.5[) then remain true for any polynomials p and q of degree 
at most n — 1 . 

We also remark that in addition to (|2.4p - (|2.5p . there are many other formulas of reproducing 
type for the kernels Kij, i,j — 1,2. We prefer to state them not all separately but will sometimes 
encounter them in the proofs below. 

2.2 Proofs of Theorems Q and [TT51 

In this section we prove Theorems 11.21 and 11.31 The proofs will use the Riemann-Hilbert charac- 
terization of Kuijlaars-McLaughlin [32] . 

We recall our assumption that the potential W is a polynomial of even degree: 

W(y) = c q y q + 0{y q - 2 ), c q > 0, q := deg W £ 2N. (2.6) 

Following [35], we define the weight functions 

/>oo 

w {x):= yi e -V(x)-W(y) + TXy dy ^ j = 0) . . . , g - 2. (2.7) 



We stack these weight functions together in the row vector 

r(x) := (w (x) ... w g - 2 (x)) . (2. 



w 



The Riemann-Hilbert problem (RH problem) of Kuijlaars-McLaughlin [32] is as follows. See 
also [TTJ [T3J [53J [25] for alternative Riemann-Hilbert problems for the two-matrix model. 

RH problem 2.2. We look for a matrix-valued function Y(z) of size q by q such that 

(1) Y(z) is analytic in C \ R; 

(2) For x G E, it holds that 

Y + (x)=Y_(x)(l ^^j, (2.9) 

where I q -\ denotes the identity matrix of size q — 1, where w(x) is defined in (|2.8p . and 
where the notation Y+(x), Y_(x) denotes the limit ofY(z) with z approaching x € R from 
the upper or lower half plane in C, respectively; 

(3) As z — > oo, we have that 

Y(z) = (I q + 0(1/ z)) diag(z", z~ n \ z' n ^), (2.10) 
where n k := [ n+q Zl~ k } , k = 0, . . . ,q - 2. 
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Similarly, we may consider the following dual RH problem. 
RH problem 2.3. We look for a matrix-valued function X(z) of size q by q such that 

(1) X(z) is analytic in C \ R; 

(2) For holds that 

X + (x)= JT_(*)(V W 7 X) )' ( 2 - U ) 

where the superscript T denotes the transpose; 

(3) As z — > oo, we have that 

X(z) = (I q + O(l/z))dmg(z n0 ,...,z n ^,z- n ). (2.12) 

Partition the matrix Y = Y(z) as 

1 q-1 

Y(z) = 1 fY ltl (z) Y h2 {z)\ (2.13) 



q-1 \Y 2 ,i(z) Y 2>2 {z)) ' 
where the partition is such that Yi i has size lxl. Similarly, partition the matrix X — X(z) as 

q-1 1 

X(z)= q-1 (X ltl (z) X 1>2 (z)\ (2.14) 
1 \X 2A {z) X 2 , 2 (z)J' 

where the partition is such that X\,\ has size (q — 1) X (q — 1). 

The solutions to the RH problems 12.21 and 12.31 are related through the formula 



X = ( XlA Xl 'A = ( Y2 ' 2 - ^ 2 >A t (2.15) 
\^2,i X 2<2 J \—Yx t2 Yi i J ' 

where the superscript ~ T denotes the inverse transpose. This was first shown by Van Assche et 
al. in @D], see also HHB]. 

Lemma 2.4. With the partitions (|2 . 13|) and (I2.14p as described above, we have that 

Y 1A (z)= Pn (z), (2.16) 
which is the biorthogonal polynomial in (jl.4l) . and 

X 2<2 {z) = -±—Q n _i( z ), (2.17) 

"•n-l 

where h n -\ and Q n -i o,re as defined in (|1.7p and (|1.13p respectively. 

Proof. The statement (|2.16[) was established by Kuijlaars-McLaughlin [32] . 

Next we prove (|2.17p . To this end we recall some of the identities of Van Assche et al. [3D] : 
it is observed in the latter paper that 

X 2 ,x(z) = (Ao(z) ... A q _ 2 (z)) (2.18) 
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for certain polynomials Aj(z) of degree at most rij — 1, j = 0, . . . , q — 2. These polynomials are 
such that the following 'type 1' orthogonality relations hold: 

/oo 
x J X 2 ,i(x)w T (x) dx = -27r% n _i, (2.19) 
-oo 

for j = 0, 1, . . . , n — 1, where ^>-i denotes the Kronecker delta. Finally, [JO] also shows that 

1 f°° 1 

X 2 ,2(» = -7^ / X 2tl (x)w T (x) dx. (2.20) 

Now we are going to use the special form of the weight functions Wj{x) in (12.7I) - (|2.8[) . Using 
(|2.18[) and an integration by parts argument as in |32) it is easy to see that 

X 2 ,\{x)w T (x) := A a (x)w (x) + . . . + A q - 2 (x)w q - 2 (x) 

/OO 
B n ^(y)e~ v ^ w ^ +Tx y dy, (2.21) 
- OO 

for a certain polynomial B n -\{y) of degree at most n — 1. Inserting this in fj2.19p - fl2.20p . these 
relations are transformed into 

x j B n _ 1 (y)e- v W- w M +rxy dx dy = -STtt^.n-i, (2.22) 



— coJ — OO 



for j = 0, 1, . . . , n — 1, and 



1 /■OO /■OO 1 

X 2 , 2 {z) = -—j j ——B n ^(y)e- v ^- w ^+ rx y dx dy, (2.23) 



— ooJ — oo 



respectively. 

Now (|2.22p shows that -B n _i is precisely the biorthogonal polynomial q n -i in fj 1 .4|) . up to a 
constant, i.e., 

B n -i(y) = C„_iq„_i(y), 
for certain C n -i S C. Then (|2.23[) implies in turn that 

X 2 , 2 (z) = / — q n -l{y)e- V{x) - W ^ +Txy dx dy (2.24) 

Z7TZ J-ooJ—OQ Z x 
C n -1 



2-ki 



-Qn-i{z), (2.25) 



where the second equality follows by virtue of (I1.13[) . 

Finally, we want to identify the normalization constant C n _i in (|2.25D . To this end, we 
substitute the expansion — ^— = 7 + ^ + §3 + ■ • ■ and the orthogonality relations fj 1 .4[) in (I2.24[) 
to obtain that 

X 2t2 (z) = -^h 2 n _ lZ - n + 0(z- n - x ), z^oo. 
2tti 

Comparing this with the asymptotics in (I2.12[) . we see that 

2m K ~ X ~ L 

Inserting this expression for C„_i in (|2.25[) . we obtain the desired formula (|2.17l) . This ends the 
proof of the lemma. □ 
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Proof of Theorems and \1.3\ Taking into account Lemma 12.41 Theorems 11.21 and 11.31 now 
follow directly from the results in [T3] and [3U] respectively; see also [TH] OH] ■ □ 



Remark 2.5. The kernels i^x.i in (jl.8|) and Ki t \ in (|2.1I) can be expressed in terms of the RH 
problem for Y(z) as follows: 



1 



2wi(xi — X2) 



(0 w(a: 2 )) y-Hasa)^!) 







and 



Ki t i(x,v) 



X — V 



1 



0) F" 1 («)y(a;) 







(2.26) 



(2.27) 



These formulas can be obtained from the results in [18] and [19] respectively. Note that if 
x\,X2 6l (or then we should replace Y in (|2.26[) (or (|2.27[) respectively) by one of its 

boundary values Y + or y_; both boundary values lead to the same formulas. Also note that, by 
means of (|2.27p (or (|2.1[1 ) and the results in [T5], it is possible to obtain determinantal formulas 
for pi 1 ' *® for any I and K; however we will not follow this route here. 



2.3 Proof of Theorem [T31 

In this section we prove Theorem ll.5l The proof will be virtually the same as the one of Akemann 
ct al. [3]; we include it for convenience of the reader. The proof will require two formulas of 
Christoffel type, see [3] © [37]. 



Proposition 2.6. Define the weight function 

w(x,y) : 



Y[(x - Xi) e" 

i=i J 



V(x)-W(y)+Txy 



(2.28) 



with I > 0. The monic biorthogonal polynomial A n {x) of degree n with respect to this weight 
function, which is defined by the biorthogonality relations 



OO /-OO 



A n (x)q(y)w(x, y) dx dy = 



(2.29) 



—00^—00 



for any polynomial q(y) of degree at most n — 1, is given by 

(p n {x\) ... p n+ l{xxf\ 

det 

1 



A n (x) 



Pn(Xl) 
\Pn(x) 



p n +l(xi) 
Pn+l{x) J 



det : 

\Pn{xi) 




(2.30) 
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Proof. Clearly A n (x) is a monic polynomial in x of degree n. So it suffices to check the biorthog- 
onality relations (|2.29[) . Observe that in the integrand A n (x)q(y)w(x, y), the factor Yii=i( x ~ x i) 
of w(x,y) in (|2.28l) cancels with the prefactor of A n (x) in (|2.30p . By the linearity of determi- 
nants, the double integral in (|2.29p can be performed cntrywise inside the last row of the matrix 
in the numerator of (|2.30[) . Then it suffices to prove that 



oo i>oo 



p n+i (x)q(y)w(x, y) dx dy = 0, i = 0, . . . , I, 



ooJ — CO 



which is a direct consequence of the biorthogonality relations. □ 
Proposition 2.7. Define the weight function 

w(x,y):= (n^-^II^-^) ] e- v ^- w ^ + ™y, (2.31) 

V=l 3=1 / 

where I > J > 0. The monic biorthogonal polynomial B n (y) of degree n with respect to this 
weight function, which is defined by the biorthogonality relations 



oo poo 



p{x)B n {y)w{x, y) dx dy = (2.32) 



— oo J — OO 



for any polynomial p of degree at most n — 1, is given by 
B n {y) := , kl+J 



det 



n j= i(y-%) 

(K\,2{x\,y\) ... Ki t2 {xi,yj) K ia {x\,y) p n +j+i(xi) ... p n +i-i(xi)\ 
XKx^jx^yi) ... K 1 , 2 (x I ,yj) K lj2 {xi,y) p n +j+i(xi) ... Pn+i-ijxj)/ 

X ~ ~ 

(Ki.iixx^yx) ... K ia (x!,yj) p n+J (xx) p n+J (xx) ... p„ + /_i(a;i) N 

det 

\Ki t 2 {xii yi) ... Ki a (xi,y,j) p n+J (xi) p n +.j{xi) ... p n+ i-i{xi)/ 

(2.33) 

Here we define Ki_ 2 (x,y) as in (|1.19p but with n replaced by n + p for an arbitrary but fixed 
number p G {J + 1,...,/}. 

Proof. Clearly B n (y) is a monic polynomial in y of degree n. Next we check the biorthogonality 
relations (|2.32l) . To this end we will take the index p in the statement of the lemma equal to 
I. Observe that in the integrand p(x)B n (y)w(x,y), the factor Oj=i (u ~ Vj) °f ^O^y) m (|2.31j) 
cancels with the prefactor of B n (y) in (|2.33p . By linearity we can then take the double integral 
in (I2.32p entry wise inside the (J + l)th column of the matrix in the numerator of (|2.33p . Hence 
the integral (|2.32p can be written as a linear combination of terms of the form 

f>oo poo 1 



/OO pQO 1 
/ p(x) ]J(x - x l )K 1 . 2 (x k ,y)e~ v ^- w ^ +Tx y dx dy 
-ooJ — oo i—1 

/oo J 
p(x) Y[(x - x t )K hl (x k ,x) dx, (2.34) 
-oo ■ 1 
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for k — 1, . . . ,1, where we used (II. 8[) and (|1.19p . But the integral (|2.34p is zero for any k = 
because of the reproducing property (I2.4[) for Kx,i. (Recall that we are taking the index 
p in the statement of the lemma equal to I.) □ 

Theorem [T3] now follows from Theorem 11.21 and Propositions 12 .61 and 12 . 7l bv an easy induction 
argument, see e.g. [6]. □ 



2.4 Properties of the kernels K^j 

Our next goal is to prove Theorem 11.61 In the present section we first collect some preliminary 
results on the kernels JQj which will be needed in the proof. 

First we discuss the kernel Ki^ 2 (x, y) = Kx i2 (x, y). Observe that by the definition (|1.19[) . the 
kernel K\ t2 (x, y) is a bivariate polynomial in x and y, of the form 

71-1 

K 1<2 (x,y) = J2 (2.35) 

i,j=0 

for suitable coefficients Ci_j G C. 

Lemma 2.8. (Reproducing property of Ki j2 :) The kernel K\ 2 is reproducing in the sense that 

Ki, 2 (x, y)p(Oe- v <0- w (">+ r *» d£ dy = p{x), (2.36) 

K x , 2 {x, y)q(ri)e- v W- w M+ Txr > dx d V = q(y), (2.37) 



OO fOO 



— ooJ — OO 



OO /■ OO 



ooJ — OO 



for any polynomials p and q of degree at most n — 1. Moreover, the reproducing property (I2.36[) 
(or (I2.37[) ) uniquely characterizes Kx )2 (x,y) over all bivariate polynomials of the form (|2.35D . 

Proof. Taking into account (|1.8[) . (jl.lip and f|l . 19[) . the reproducing properties (|2.36p - (|2.37[) are 
just a restatement of (|2.4p and (|2.5I) respectively The claim about uniqueness is also easily 
shown, see also [19]. □ 

Next we discuss the kernels and K 22 . We will need the following result. 

Lemma 2.9. (Vanishing properties of and K 22 :) The kernel -Ki,i satisfies the 'vanishing 
property' 

OO />OG 

Ki,i (£, v)q{r 1 )e- v ^' w ^ +T ^ d£ drj = (2.38) 



OO J — OO 



for any v G C \ IR and all polynomials q of degree at most n — 1. Similarly, the kernel i^2,2 
satisfies the vanishing property 



OO />OG 



— ocJ — OO 



^ 2 , 2 K 77)p(£)e- F ^- w W> +TS " ^ d7] = o, (2.39) 
/or any id£C\1 and aZ/ polynomials p of degree at most n — 1 . 
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Proof. This lemma can be obtained from (12.27)) and the vanishing property in |19) . For com- 
pleteness, we include a direct proof. Using (|2 . 1 [) we have 



r°° r°° ~ r°° r°° n ( n\p~ v (£)- w ( r i)+ T t r i 
[ K ltl ^v)q( V )e- v ^- w ^^d V = - [ I q{v)e 

K 1 ^x)q{i 1 )e^ v ^- W ^ +T ^ 



-OC«/ — OO J — OOJ — OO 

roo poo 



OO J — OO" 



Now we use the following formula of reproducing type: 

y* OO 

K hl {^x)q{r ] )e- V ^ w{rl)+T ^ d£ dr, = / g(r ? ) e -^)- w W+™7 dn, (2.41) 



OO f OO 



— OO tJ — OO 



for any polynomial q of degree at most n — 1. (As usual, by linearity it suffices to check (|2.41[) 
when a equals one of the biorthogonal polynomials qt for i = 0, . . . , n — 1, in which case it follows 
easily by the biorthogonality relations.) Using (|2.41[) in the term with the triple integral in 
(|2.40[) . we see that both terms in the right hand side of (|2.40[) cancel each other and so (12 .38[) 
follows. The proof of (|2.39p is similar. □ 

Lemma 2.10. (Asymptotic behavior of biorthogonal polynomials and kernels:) We have 

p n (x) = x n + 0(x n - 1 ), x^oo, (2.42) 

q n (x) = y n + 0(y n ~ 1 ), y -> oo, (2.43) 

P n ( w ) = hlw' 71 - 1 + 0(w- n - 2 ), w^oo, (2.44) 

Q n (v) = hlv-* 1 - 1 + 0(v- n - 2 ), v^oo, (2.45) 

and 

K hl (x,v) = h^Q^vy 1 - 1 +0(x n - 2 ), x-^oc, (2.46) 

Ki,i(x,v) = -pnix^- 71 - 1 +0(v- n - 2 ) , v^oo, (2.47) 

Ki, 2 (x,y) = ft-fi« n -i(l/)i n ~ 1 +0(x n - 2 ), x^oo, (2.48) 

K h2 (x,y) = /i i ;! 1 p„_ 1 (x)y"- 1 +0(y"- 2 ), y -> oo, (2.49) 

and 

# 2 ,iK«) = -^("K^+O^ 2 ), w^oo, (2.50) 

#2,i(u>,«) = -PnH^'+Off^ 2 ), v^oo, (2.51) 

K 2t 2(w,y) = -q n {x)w- n - 1 + O (w~ n - 2 ) , w -)• oo, (2.52) 

tf 2)2 («M/) = d^-iHr'+Of!/" 2 ), 2/^oo- (2-53) 

Proof. Equations (|2.42p - (|2.43D are obvious by definition. Let us check (|2.44D . By substituting 

the series = V" 00 -, ^— ?- in (11.121) we find that 

00 1 / n OO /» OO 

Now by the biorthogonality relations we have 

f OO pOO 



OO J — OO 
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for any i = 1, . . . ,n + 1. So we obtain (I2.44j) . In a similar way one checks (|2.45l) . 

Equation (g^jg]) follows by virtue of (fl~Tg)) and (PQ2")) . Next we check (j!T4T)) . By substituting 
the expansion ^| = tT" m (Ell) we obtain 

~ 1 °° 1 / f°° \ 

Jif 1)1 (x, w ) = — / (x-Or^M^O^ ■ (2.54) 

X V i=i V \J—oo / 

Now 

/oo 
(x-Of-^i.iM d£= [(a?-Or _1 ] €=ie = 0, ifi = l J ...,n-l, 
-00 

because of the reproducing property (|2.4I) . For i — n we have 

3 /.OO / 1 \ 



f°° 1 

J —00 "'n 



The expression between brackets in the first term in the right hand side of (I2.55[) is nothing but 
the kernel ^1,1(2:, £) with n replaced by n + 1 in fj 1 .8[) . so this integral is zero as before. We then 
have 

o ' ' n J — 00 

1 />00 />00 

= -TiPnix) / (z - C)r" 1 9n(r?)e-^ ) - w( ' 7)+T ^ de A? 

n •/ — 00J — 00 

because of biorthogonality. Using this relation in (|2.54l) . we obtain (|2.47p . 

Equations (|2^8]) - ([2~4"5|) are immediate by the definition (fTTT^j) . Equations ([^32]l - (|2T55|l can 
be obtained similarly as before. 

Finally, let us check (|2.50[) . By substituting the expansion — \p~ m (HUHl) we 

obtain 

K 2 , 1 (w,v)=Y / — —A rf- 1 K %1 {r,,i)di 1 )d^. (2.56) 

i_i ^ J —oo " S \J — oo / 

By (fTTU)) we see that 

/OO /> OO />OC 1 

-00 J —00 J— 00 i=o j 

(2-57) 

Now we claim that 

/oo 1 /*oo 

?(»?) E TS P i{l)Qi@ *? = / g^e"^)-^^" dry (2.58) 
-co j =0 n j J -00 

for any polynomial q of degree at most n — 1. (Once again, by linearity it suffices to check 
(|2.58p when g equals one of the biorthogonal polynomials qi for i = 0, . . . , n — 1, in which case it 
follows directly from the biorthogonality relations.) Using (12.58)) with q{rf] = rf^ 1 , we see that 
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if i = 1, . . . , n then both terms in the right hand side of (|2.57p cancel each other. For i = n + 1 
we get 

/oo re — 1 

(/■oo n , \ „oo i 

J-oo j=Q 11 j J J-oo n n 

= ( r n n e- v ^- w ^ + ^ d V ) - r r v n Pn ( X ) e - v ^- w ^+^ ^ ^ 

\J Co / lf, n J — oo*/ — oo 



where we used (|2.58l) (with n replaced by n + 1) and the biorthogonality relations. Substituting 
this in (j2~5Tf we find that 



/OO 
f? i-1 -K2,l(r?,£) *7 = -^,n+lQn(6) 
-oo 



for any i = 1, . . . , n + 1. Using this in (12. 56ft we obtain 

#»,!(«>> «) = -W-"- 1 r dC + 0( W - n - 2 ) = -W-^Qniv) + 0( W ~ n ~ 2 ), 

J-oo V-£ 

which is (|2.50p . Equation (I2.51[) can be obtained similarly. □ 

Finally, we will also need the following additional properties of the kernels Kij. 
Lemma 2.11. (Integral formula relations between the kernels Ki^:) We have 

/OO ^00 i 
/ Pn (0e- ni) - mr,)+T$v d£dn, (2.59) 

K 2 ,2{w,y) = — I" r K 1 . 2 (^y)e~ m) ' WM+T ^ l]i ^ 1 - d£ d V , (2.60) 

y-w j_oo J-oo w-v 

and 

/OO /"OO 1 
/ K 1 , 1 (Z,v)e- v ^- w( ^+^ d£dr,. (2.61) 

Proof. Equation (12.591) is just (|1.12l) . Equation (|2.60p follows from the definitions in (|1.17l) . 
([TTTTj) and (TTTO1) . Finally, equation (|2lflT) follows by (TOO)) . (fTTTHjl and (|TT2"|) . □ 

2.5 Proof of Theorem 11.61 



In this section we prove Theorem 11.61 To this end we use Theorem 11.51 together with a mech- 
anism to transform an external source in the matrix M\ (or M%) in the numerator, into an 
external source in the matrix Mi (or M\ respectively) in the denominator. That is, we will 
apply transformations of the form 

(I,J,K,L) ^ (I -l,J,K,L + l), (2.62) 
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or 

(/, J, K, L) (I,J — 1,K + 1, L). (2.63) 

It will then suffice to show that each of the objects in (|1.12[) - (I1.17I) transforms in an appropriate 
way under these transformations; this will be achieved by virtue of Lemmas I2.8H2.1T1 

As in [SJ [TO], we will also need an appropriate use of partial fraction decomposition. For 
given complex numbers yi, . . . , yj, Wi, . . . , wl+i £ C, the partial fraction decomposition which 
is of interest to us is of the form 

n ^ l( V j) = - £5±? - + Y — + P(n). (2-64) 

Here c%,. . . ,cl+i G C and P is a polynomial of degree J — L — 1. (We put P = when 
J-L-KO.) 

Note that (|2.64l) also implies the more complicated partial fraction decomposition 

UUv - w) _ o L+1 (v - y 3 ) , £ + m _ p(pj)) 



(ry-ui L+ i)nf = i(?7- w z ) K+i -Vj)(v- w L+i) (™l ~ Vj)(v ~ ™l) 



(2.65) 



for any fixed j G {1, . . . , J}. We leave it to the reader to obtain (|2.65[) from (|2.64j) . 
Now we are ready for the proof of Theorem 11.61 

Proof of Theorem \1.6{ By symmetry we can assume without loss of generality that I — K > 
J — L. We will show that the formula in (|1.28[) is compatible with any transformation of the 
form (gllg) . 

Consider the weight function 

W( x ,y) = ~ Xl) ^ liV - V A e -V{*)-W{ V )+T X y_ (2 g6) 

Y[ k =i( x - v k)Y[i=iiy - wi) 

(Note that there are only I —1 factors X{.) Denote by A n (x) the monic biorthogonal polynomial 
of degree n with respect to this weight function, defined by the orthogonality relations 



oo /-oo 



A n (x)q{y)w(x,y) dx dy = (2.67) 



oo./ — oo 



for all polynomials q of degree at most n — 1. Theorem 11.21 implies that 

p [i,j,k,l}^ x ^ ,xi;y!, . . . , yj\ v%, . . . ,v K ;wx, . . . , w L ) = A n (xi), 
while from Theorem II .31 and (|1.12D it follows that 

p[I-l,J,K,L+t\( Xi ^ _ _ . )2;j _ i; y lj ...,yj; Vl ,.. .,V K ;W!, . . .,W L+1 ) 



X ' V w(xj,Tl) 

A n -i(xi) dxi drj, 



hi, J -ooJ-oo " ^ "w L+ i-rj 



where h„_i is a constant which depends on each of the numbers Xi,yj, Vk, wi, but not on Wl+i- 
By combining these two formulas, one gets 

p[I-l,J,K,L+l)( Xi ^ _ _ , iXl _ i;yi} ...,yj; Vl ,.. .,V K ;Wi, . . .,W L+1 ) 

1 />00 r CO ~ / \ 

P' n _ 1 '{xi,...,x I ;yi,...,yj;vi,...,v K ;wi,...,w L )— - dxi dt]. 

(2.68) 



h 2 n _ x J-ooJ-oo WL + 1 - V 
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Now assume, by induction, that (ll.28|) holds for Pn' J ' K ' L ^. By substituting this expression (with 
n — 1 instead of n) for p^^ K ' L ^ i n the right hand side of (I2.68[) . we get 

p\f-\,J,K,L+l]^ _ _ _ )X/ _ i;2/l; ...^jjui,...,^^!,.. ■ ,WL+l) 



Co 



x det 



/ ^1,1(0;!, vi) 

Ki, 2 {xi,yj) 

Pn-l+J-L{x\) 
\Pn-2+I~K(Xl) 



Ki,i(xi,v K ) K2,i(wi,v K ) 

Ki t2 (xi,yi) ^2,2(^1,2/1) 

Ki,2(xi,yj) K 2 , 2 {w\,yj) 

Pn-1+J-L(xi) Pn-1 + J-L(wi) 



K 2jl (w L ,Vi) \ 

K2,i(wl,v k ) 
K 2 , 2 (w L ,yi) 

K 2 , 2 {w L ,yj) 

Pn-1 + J-L(W L ) 



dxi drj, 



Pn- 2 +I-K{xi) P n -2+I-K{wi) ... Pn- 2 +I-K{w L )/ 



(2.69) 



where Co is a new constant which depends on each of the numbers Xi, yj, v k ,wi, but not on wl+x- 

pr/-l/ _ -i _____ 

We will work this out. Note that the factor ffif 1 , s in the definition of w; (27,77) in (I2.66[) 

1 lk=l \ x I~ v k ) 

cancels with the prefactor in the integrand of (|2.69p . By linearity, we can then apply the double 
integration entry wise inside the 7th column of the matrix in (j2.69[) . Then the entries in the 7th 
column of (12.69)) transform into expressions of the form 



00 poo 



Ki A (xi,v k ) 



— OCJ — OO 



—i < 

Ui=i(v-wi) wl+1 -V 



for k = 1, . . . , K, 



00 poo 



oo<J —00 



Ki 2 (xi, yj) L e 



-V(,xi)-Wb)+rxiV dxi drj) (2.70) 



vixri-ww+Txtr, dxi dr?) (2.71) 



for j = 1, . . . , J, and 



.... v. --00 Ui=i( 7 l- wi) -V 



(2.72) 



for p = J — L, . . . , I — K — 1. 

Next, we substitute the partial fraction decomposition (|2.64p in (|2.70p and (|2.T2[) . and we 
substitute (|2.65[) in (12 . 71|) . Doing this for each of the entries in the 7th column of (12 .69|> . the 
determinant can be split in a sum of three terms, 



B x +D 2 + D 3 , 

corresponding to the three terms on the right hand sides of (|2.64|l — (|2.65p . 



(2.73) 
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For the first term D\ (which is obtained by selecting the terms CL+1 and -, CL+1 ^ 

— , ... 77— tuz,+i {w L + l~Vj){v-w L + 1 ) 

in (|2.64p and (|2.65|) respectively), the expressions (|2.70[) - (I2.72[) transform into 



i,i(a:/,t; fc ) e -v(^)-w(r,)+r Xlv dxj ^ 



'■L+l I K 

-ooJ-oo W L+1 -T] 



(2.74) 



for k = 1, . . . ,K, 



V-Vj 



roo poo 

>-L-i j / Kx, 2 (xi,yj)- — 

t-ooJ-00 (w L+1 - yj)(w L+ i - 77) 



for j = 1, . . . , J, and 



CL+l 



00 POO 



Pn-X+p{xi) e ~V( Xl )-W M+ r XlV ^ ^ 



— OO </ —OO 



WL+1 - V 



(2.76) 



for p = J - L, . . . , I - K - 1. Comparing this with Lemma |2~TTI we see that ([2773]> - (p7TB]) 
are nothing but cl+i times K 2 ,i(wL+i,v k ), K 2 ^(w L+ i,yj) and P„-i+ p (wl+i) respectively. The 
factor cl+i can be taken out of the 7th column and be put in front of the determinant. Note 
that 

Uj=i( w L+i ~Vj) 

CL+1 = W~i \' 

a fact which is easily checked from (|2.64p . 

For the second term D 2 in (12.73)) (which is obtained by selecting the terms J2i=i r)-w l anc ^ 
Yd=x ( w ^-y^l Wl ) in (|^M]) -I[2~55 | respectively), we obtain similarly that (|277D)l - (|2~?2")) trans- 
form into a linear combination oi K 2r i(wi,Vk), K 2 , 2 (wi,yj) and P n -i+ p (wi), I = l,...,L. But 
then the 7th column in (|2.69[) is a linear combination of columns I + 1 , . . . , I + L and so the 
determinant D 2 vanishes. 

Finally, the third determinant -D3 (corresponding to the polynomial parts P(rj) and P{rf) — 
P{yj) in (|2.64p ~ (|2.65p respectively) vanishes as well, by virtue of Lemmas 12.81 and [231 

Summarizing, we showed that (|2.69p equals 

D \i-xjKL+X[ t \ n n/=i(^i+i -Vj) 
Pk ' i {x 1 ,...,x I - 1 ;y 1 ,...,yj;v 1 ,...,v K ;w 1 ,...,w L+1 ) = C\ " 

1L=i( w £+i - w 



x det 



Ki,i(xi,Vk) ■•■ Ki,i(xi-i,vk) K2,i(wi,v k ) 

-^1,2(21,2/1) ■•■ K 2 ,i(xi-i,yi) K 2<2 (wi,yi) 

Kia(xi,yj) ... K 2A (x I ^ 1 ,y J ) K 2 , 2 {w x ,yj) 

p n -l+J-L(xi) ... p„-l + J- L (xi-l) Pn-i + j-LiWx) 



K 2 , 2 (w L +i,yi) 
K 2 , 2 {w L+1 ,yj) 

Pn-l + J-L(w L + l) 



\Pn- 2 +I-K(xi) ... Pn- 2 +I-K{xi-l) P n - 2+ I- K {Wl) ■ 

where C\ — ±Cq depends on each of the numbers Xi, yj, v k , toj, but not on wl+i 



P n - 2+ I-K{wL+l)J 

(2.77) 
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To find Ci, we compute the leading order behavior of (|2.77[) in Wl+V, this can be done by 
replacing each of the entries in the last column of (|2.77l) by their dominant term in wl+i as in 
Lemma \2. 101 (More precisely, if/ — K > J — L then only the entry P n -i+j-i(wL+i) in (12.77)) 
contributes to the dominant term, while i£I — K=J — L then there is a contribution from all 
entries in the last column of (|2.77p .) Taking into account that 



J L+1' L n 



,yj;v!, ...,vk;wx,.. .,w l+1 ) 

, . , xi-i; 2/1, . . . , 2/j; wi, • • • , v K ; Wi, . . . , w L ) + O (w^ 1 ) , 



as wl+i — > co, a fact which trivially follows from (| 1 . 24[) . we then obtain an expression for C\. 
In a similar way we can compute the leading order behavior of (|1.28[) in xi, i.e., we can write 



P 



[I,J,K,L] 



in the form 



a 



Uk=i( x i 



Wfc) 



K 1:2 {xi,yi) 



x dct 



Ki l2 (xi,yj) 

Pn+J-L(Xl) 
XPn+I-K-l(xi) 



Ki,i(xi,vk) 

K 2 ,i(xi,yj) 
Pu+j-l(xi) 



Kj.i(wi,v k ) 
K 2 .2{w 1 ,y 1 ) 

Pn+,J-L{wi) 



Pn+I-K-l(xi) P n+ I- K -l(wi) 



K 2i i(w L ,vi) \ 

K 2:2 (w L ,y 1 ) 

K 2 , 2 (w L ,yj) 
Pti+j-l(wl) 

Pn+I-K-l(w L )J 



, (2.78) 



and then we can replace each of the entries in the 7th column of (|2.78D by their dominant term 
in xi as in Lemma 12.101 (More precisely, if J — K > J — L then only the entry p n +i-K-i(xi) 
in (|2.78l) contributes to the dominant term, while if/ — K = J — L then there is a contribution 
from all entries in the Jth column of (j2.78|l .) Taking into account that 



P. 



I.J.K.L 



] (x 1 ,...,x I ;y 1 ,...,y J ;v 1 ,...,v K ;w 1 ,...,w L ) 



= x nplX-l,J,K,L] ixu ... tXl _ v 



2/1, ■ 



, yj; vi, . . . , v K ; wi, . . . , w L ) + O (x" x ) , 



as xi —¥ co, we then obtain an expression for C 2 in (|2.78l) . A straightforward calculation now 
shows that the so obtained expressions for C\ and C 2 in (|2.77[) and (I2.78|) are equal to each other, 
up to a shift n i Y ti 1 of the index, and a factor (— 1) I+J+K+1 /hf l+J _ L _ 1 . So the constants C\ 
and C 2 are related precisely in the way that is required for applying the transformation (|2.62[) 
to the prefactor in (jl.28|) . 

Summarizing, we have shown that (|1.28[) is compatible with any transformation of the form 
(12.62)) . In a similar way one shows this for transformations of the form (12.63)) . Then applying the 
transformations (|2.62[) - (|2.63j) repeatedly, and using Theorem 11.51 as induction basis, we obtain 
(OS)) in its full generality. □ 



Remark 2.12. An alternative proof of Theorem 11.61 can be obtained by establishing formulas of 
Christoffel-Uvarov type [H [35] ■ To this end one can use similar ideas as in the proof above. 
We preferred the given proof since it is probably shorter. 
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3 Applications 

In this section we give two applications of our results. Here we follow ideas for the 1-matrix 
model. 

3.1 Generating function for averages of products of traces 

In this section we show how the results in this paper allow to obtain expectation values of 
products of traces of the form 



or equivalently 



V J f n^Mf)}]^) e Tr(-V(M 1 )-W{M 2 )+rM 1 M 2 ) ^ m ^ 



i=l j=l 



r •••/ ntec nte^ 

n n 

II (e-^)e-^)) A(A)A( M ) det(e^ )& =1 J] (rfA, d W ) , (3.1) 



x 

i=l 



with exponents mi, . . . , m/;rii, . . . , nj e NU{0}, see (| 1 . 3|) . Here we follow Bergere [5]. We start 
from (TQ5)) with I = K and J = L: 

i r r -q - w ^ ( nLifa - **) 



n n 

x 1] (e- v ^e~ w ^) A(A)A( M ) det( e - A ^ )™ - =1 J] (d\ t dm) . (3.2) 



Applying the operator 



to (13. 2p leads to the (multi-variate) Cauchy transform 



n \ J / n 



Znl-co'" Loo n (£ ,-, aJ n (E , /; - Mfc 



II (e-^^e-^)) A(A)A( M ) det(e^ )? . =1 JT (dA, d W ) , (3.4) 



x 

i=i i=i 



where we assume that all Xi,yj 6 C \ R. Note that the large Xi, yj expansion of (|3.4[) is a formal 
power series whose coefficients are the expectation values of products of traces (|3.1j) . 
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On the other hand, Theorem 11.61 with I = K and J = L shows that (I3.2[) equals 

nlfe=i i x i - v k) Uj,i=i(yj ~ m) 

y l )(w i - Wj) 



Ui< l<3 <i( x 3 - %i)( v i - v j) Ui<i<j<j(yj - Vi)( w i - w j) 



x dct 



K ltl (x!,vi) ■■■ Ki t i(xi,vi) Ki,2{xi,yi) ... K 1)2 (xi,yj) 
K 2 ,i(wi,vi) ... K 2 ,i(wi 1 v I ) K 2>2 (wi,yi) ... K 2t2 (wi,yj) 



(3.5) 



\K 2 ,i{w J ,v 1 ) ... K 2 .i(w Jl v I ) K 2<2 {wj,yi) ... K 2<2 {wj,yj)j 
Applying the operator Q3.3[) to this determinant, and making a small calculation, one gets 



det 



( K\ t \{x\,x\) 
K 2> i(yi,xx) 



Ki t i(xi,xi) Ki t2 (xi,yi) 

K lt i(xi,xi) Ki,2(xi,yi) 
K 2t i(yi,xi) K 2>2 (yi,yi) 



\K 2A (y J ,x 1 ) ... K 2 ,i(yj,xi) K 2>2 (yj,yi) 
where the entries on the main diagonal are defined by 



Ki, 2 (x!,yj)\ 

K l!2 (xi,yj) 
K 2)2 (yx,yj) 

K 2}2 {yj,yj)J 



(3.6) 



and 



n— 1 ^ 

Ki,i(x,x) = ^ T2Pi( x )Qi( x )> 
t=0 i 



n — 1 ^ 

K 2 . 2 (y,y) = ^2 ToPi{y)qi{y), 



respectively; compare with (|1.18[) and (11 .2 1|) . So f|3.6[) yields a determinantal formula for the 
multi-variate Cauchy transform Q3.4p . By the above discussion, this is also the generating function 
for averages of products of traces (|3.1[) . 

3.2 Eynard-Mehta theorem for correlation functions 

In this section we show how Theorem 1 1 . 6 1 can be used to give an alternative proof of the Eynard- 
Mehta theorem for correlation functions in the two-matrix model. Alternative proofs of this 
formula, which are applicable for the more general model of random matrices coupled in a chain, 
can be found in [HI [10 E3 E3 [38] . 

The correlation function Ri t j is defined from (jl.3l) by 



Rj,j{Mi ■ ■ ■ t Aj; Hi, . . . , — 



(n-I)\ (n- J)! 



1 

x — 

Z. 



n *J — oo 



n 



<r v ^e- w ^) A(A)A(/i)dfit(e TA *W)"j=i II dX * II ( 3 ' 7 ) 

»=I+1 j = J+l 
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The Eynard-Mehta theorem |24| then asserts that 



-R/,j(Ai, . . . , \i\ Hi, ■ . ■ , hj) 
( ^i,i(Ai, Ai) 



det 



^i,i(A/,A!) 
-^2,1(^1, Ai) 

\k 2 ,i(hj, Ai) 



Ki,i(Ai,Ai) Xi, 2 (Ai,mi) 

Kia(Xi, A/) K lj2 {Xi,fii) 

K 2 a(hi,\i) ^2,2(^1,^1) 

K 2 ,i{lJ.ji A/) K2.2i.lJ.. j ^1) 



K lt2 (Xi,fJ,j)\ 
Ki. 2 (Xi,Hj) 

K 2: 2(fJ.J,Hj)J 



(3.8) 



We now establish this formula by using the formulas in Section 13.11 Here we follow [TB] . We 
start from the multi-variate Cauchy transform in (|3.4[) . Applying the operators 



f !->■ lim ( f L =A ■- 



f\xi— Xi+ei) 



subsequently for i = 1, . . . , I and then 



1 



for j = 1, . . . , J to (|3.4[) . we obtain by the Stieltjes-Perron inversion principle (37) precisely the 
correlation function (I3.7[) . On the other hand, applying these same operations to (|3.6p and 
using again the Stieltjes-Perron inversion principle leads to the right hand side of (|3.8p . This 
establishes (I3.8[) . 



4 Concluding remarks 

The two-matrix model is a particular instance of a more general model, sometimes referred to 
as random matrices coupled in a chain, see e.g. [21 1151 1241 1271 1331 1351 138] . The Eynard-Mehta 
theorem for correlation functions can be formulated for this more general model. 

In view of this observation, it is natural to ask whether the results in this paper can be 
extended to the more general model of random matrices coupled in a chain. That is, one may ask 
whether the averages of products and ratios of characteristic polynomials in this model can still 
be written as determinants built out of (transformed) Eynard-Mehta kernels and biorthogonal 
polynomials. A little thought reveals that such a result, if it exists, should be a non-trivial 
extension of Theorem 11.61 except maybe for special configurations of the external sources. This 
is an open problem. 

Another question of interest is whether the results in this paper have an analogue for the 
Cauchy two-matrix model in |12) . 
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